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Deriving the 3-dimensional volume density distribution from a 2-dimensional light 
distribution of a system yields generally non-unique results. The case for nearby sys- 
tems is studied, taking into account the extra constraints from the perspective effect. It 
is shown analytically that a new form of non- uniqueness exists. The Phantom Spheroid 
(PS) for a nearby system preserves the intrinsic mirror symmetry and projected asym- 
metry of the system while changing the shape and the major-axis orientation of the 
system. A family of analytical models are given as functions of the distance (Do) to 
the object and the amount (7) of the PS density superimposed. The range of the ma- 
jor axis angles is constrained analytically by requiring a positive density everywhere. 
These models suggest that observations other than surface brightness maps are re- 
quired to lift the degeneracy in the major axis angle and axis ratio of the central bar 
of the Milky Way. 
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1 INTRODUCTION 

Deprojection of galaxies from the observed light distribution on the sky plane to the intrinsic 3-dimensional volume luminosity 
distribution is one of the basic problems of astronomy. It is common knowledge that the deprojected results are generally 
non- unique because of the freedom of distributing stars along any line of sight. The best example is that a round distribution 
in projection may correspond to any intrinsically prolate object pointing towards us. Likewise we cannot tell from observed 
elliptical isophotes whether the object is an intrinsically oblate bulge or a triaxial bar if both are edge-on and at infinity 
(Contopoulos 1956). 

The perspective effect of nearby triaxial objects, however, does make them appear different from an oblate object. The 
best example for this is the famous left-to-right asymmetry of the dereddened light distribution of the Milky Way bulge when 
plotted in the Galactic (l,b) coordinates, which led Blitz & Spergel (1991) to conclude that the Galactic bulge is in fact an 
almost edge-on bar, pointing at an angle from the Sun. They divide the Galaxy into a left and a right part with the I = 0° 
plane, which passes the Sun-center line and the rotation axis of the Galaxy. When folded along the / = 0° line the surface 
brightness map 1(1, b) is decomposed into two independent maps: an asymmetry map [1(1, b) ~ I(—l,b)] /2 by subtracting the 
I < 0° side from the / > 0° side, and a symmetric map [1(1, b) + I(—l, b)] /2 by adding up the two sides. The signal in the 
asymmetry map, they explain, is because the right hand side (I > 0°) of the bar is nearer to us and the perspective effect 
makes it appear slightly bigger than the left hand side (I < 0°). A simple sketch of the geometry is shown in the top diagram 
of Fig. § 

The perspective effect allows Binney & Gerhard (1996) and Binney, Gerhard & Spergel (1997) to derive a non-parametric 
volume density distribution of the inner Galaxy. The key element in their method is to impose mirror symmetry for the bulge 
part of the Galaxy to allow for a triaxial bar, and central symmetry for the disk part so to allow for the spiral arm. Although 
shifting material along the same line of sight does not alter the isophotes, it spoils the symmetry of the object. Their numerical 
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experiments suggest that the COBE map is consistent with a range of bar models with the bar major axis pointing some 
15° - 35° from the Sun. 

Unfortunately there is a lack of analytical studies of general non-uniqueness for nearby objects. This is compared to a 
series of papers exploiting the analytical properties of konuses in the deprojection of external axisymmetric systems (Rybicki 
1986, Palmer 1994, Kochanek & Rybicki 1996, van den Bosch 1997, Romanowsky & Kochanek 1997). A small amount of 
konuses, as christened by Gerhard & Binney (1996) for a well-studied class of artificial density models with zero surface 
brightness, can be added to a galaxy at infinity without perturbing its isophotes or creating negative-density zones. 

This paper is a first attempt to give an analytical description of the degeneracy in deprojecting the light distribution of 
a general nearby system. To ensure our arguments are not diluted by the necessary mathematics, we will first introduce the 
so-called phantom spheroidal (PS) density in §2 and describe its effect on non-uniqueness in words and illustrations. The more 
mathematical aspects of the problem are given in §3-§5, where we describe the properties of PS densities and how to generate 
them. We will then briefly discuss the relations to the Milky Way bar in §6, and relations with well-known non-uniqueness of 
extragalactic objects in §7. We conclude in §8. Some additional results on the major axis angle and a generalization of the 
phantom density are given in the Appendix. 

2 PHANTOM SPHEROID: A COMPROMISE BETWEEN MIRROR SYMMETRY AND 
PERSPECTIVE EFFECT 

A phantom density is a model with both positive and negative density regions and a zero net surface brightness to an observer 
at some distance away. An example of a phantom density is illustrated in the middle diagram of Fig. |l|, where we tailor 
the radial profile of a stratified ball-shaped bulge in such a way so to have the same angular size and projected intensity as 
a uniform cigar-shaped bar placed end-on. Subtracting the ball-shaped bulge from the cigar-shaped bar yields a Phantom 
Spheroid (PS). The thin dark ring in the bulge is a density peak, corresponding to the line of sight to the far edge of the cigar- 
shaped bar, a direction where the depth and the projected intensity of the bar are at maximum. The fall-off of density towards 
the edge of the bulge corresponds to the ever-decreasing depth of the cigar-shaped bar with increasing impact parameter of 
the line of sight. 

The problem with this simple cigar model is that if we superimpose this unphysical component on any physical density 
distribution, say, an ellipsoidal bar with a Gaussian radial profile, we get a new density with generally twisted density patterns. 
So while the added cigar component is invisible from the observer's perspective, it generally spoils the mirror symmetry of 
the system. 

The exception is, as shown in the lower diagram of Fig. [j], when the superimposed end-on cigar is a clone of the original 
cigar, in which case the final density of the twin should have mirror symmetry with respect to a new plane (the dotted line) 
which is just in between the Sun-center line and the major axis of the original cigar. In this case the mirror symmetry is 
preserved, only the symmetry plane is rotated. Subtracting the round bulge has no further effect on the symmetry plane, but 
will take away any trace of transformation in the map of the integrated light. Thus a PS which rotates the symmetry axes 
and is invisible in projection can be obtained by placing the original prolate cigar end-on and then subtracting off a round 
bulge. 

We can apply the same trick to oblate systems, such as disks, because a face-on disk mimics a spherical distribution 
the same way a prolate cigar bar does. So suppose the original model is a disk tilted at an angle a from the line of sight. 
Adding a face-on clone of the disk will put the mirror plane of the twin disks at an angle ^. By subtracting off a spherical 
model we can take away any effect of the face-on clone in the projected brightness map. In fact this kind of non-uniqueness 
applies any spheroidal (i.e. prolate or oblate) distribution with any radial profile since our arguments about non-uniqueness 
is independent of the aspect ratio and density profile of the bar. For example, the original model may consists of a small 
spheroidal perturbation on top of a large positive spherical component, both with general smooth radial profiles. 

Despite these generalizations, the models are restricted in the sense that the models come only as twins with identical 
projected brightness distribution. This is compared to non- uniqueness in the extragalactic systems, such as konuses, where 
we find a family of models with indistinguishable surface brightness maps with a tunable amount of konuses. The models here 
are also likely to have large zones of negative density because we subtract off a significant amount of matter with a spherical 
distribution. This problem can be softened if our original model has a large positive smooth spherical component to start 
with. 



3 METHOD FOR GENERATING GENERAL PHANTOM SPHEROIDS 

Let's first set up a rectangular coordinate system (x, y, z) centered on the object with the x axis pointing towards the observer. 
In this coordinate system the observer is at (-Do, 0,0), where Do is the observer's distance to the object center. From the 
observer's perspective the surface brightness map of the system is most conveniently specified by two angles (/', &'), which are 
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Figure 1. Top panel: the xy cross section of a uniform cigar-shaped bar at distance Do with its major axis tilted at an angle a 
counterclockwise from the Sun-object center line. The coordinates of a point in the xy plane can also be prescribed by its line of sight 
distance D, the impact parameter to, and the distance r to the object center. Middle panel: a phantom spheroid (PS). Bottom panel: a 
superposition of the above two bodies. Several line-of-sight pathes from the Sun are drawn in thin solid lines. The PS is made by putting 
a clone of the cigar-shaped bar to an end-on geometry, and subtracting from it a stratified ball-shaped bulge of the same angular size 
and projected intensity. Adding such a PS to the cigar-shaped bar in the top diagram has no effect on the latter's projected intensity, 
but rotates the mirror plane of the spatial distribution to the dotted line in the middle at an angle ^ from the Sun-center line. 



the equivalent of the Galactic coordinate system with the b' = z — plane being the xy plane and the I' = y = plane passing 
through z-axis of the system. In these coordinates the system center and anti-center are in the directions (l',b r ) = (0,0) and 
(tt,0). A point (x,y,z) in rectangular coordinates can thus be specified by the distance D to the observer along the line of 
sight (l',b') with 

x — Do — D cos b' cos l' , y = Dcosb' sin/', z = Dsinfc'. (1) 
The impact parameter w with respect to the object center is given by 

w = D ^(sin b'f + (cosfo' sinl') 2 (2) 
for the line of sight (/', b'). The distance r to the center of the object at (0, 0, 0) is given by 
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r = y/x 2 + y 2 + z 2 = y/A 2 +w 2 , A = D - yf D 2 - w 2 , (3) 

where A is the offset distance from the tangent point. The rectangular coordinates can also be expressed in terms of w and 
D with 



D 



x = Do-D x ll-^, y/y 2 + z 2 = —w. (4) 



Most generally our phantom spheroid F(x, y, z) is made by subtracting from an end-on prolate distribution P(|a;|, yf y 2 + z 2 ) 
a spherical bulge S(r) with matching surface brightness, i.e., 

o 

F{x,y,z)dD = 0, F(x, y, z) = P(\x\, ^r 2 - x 2 = yj y 2 + z 2 ) - Sir), (5) 

where f^°dD is an integration along any line-of-sight direction, say, (l',b r ). Clearly F(x,y,z) is an even function of x, y and 
z with 

F{x,y,z) = F{x,-y,z) = F(x,y,-z) = F(x,-y,-z) = F(-x,y,z), (6) 
and in fact has rotational symmetry around the :r-axis. The total luminosity of a phantom Lps is given by 

Lps= [d 3 rF(x,y,z)=L P -L s , L P = [d 3 rP(\x\, y/y 2 + z 2 ), L s = fd 3 rS(r). (7) 



As we will show later a phantom spheroid can have a net luminosity because the total luminosity of the prolate component 
Lp is only approximately that of the spherical component Ls with the difference being infinitely small when the observer is 
infinitely far away from the object. 

Let J(w) be the light intensity of the prolate distribution P(|a;|, yj 'y 2 + z 2 ) integrated over a line of sight with impact 
parameter w, including both the forward direction and the backward direction, then 

/oo 
P(\x\,y/y 2 + z 2 )dD, (8) 
- OO 

where x, y, z can be expressed in terms of w and D using eq. (j3|) and eq. (Q). Once J(w) is computed from the integration of 
P(\x\, \/y 2 + z 2 ), our task is to find a spherical bulge S(r) such that 

J(w) = / S(r)dD. (9) 



This can be done using the Abel transformation 

S (r) = -lpy f (10) 

TV J r dw y/ W 2 _ r 2 

The inversion uses effectively a variation of the well-known Eddington formula for deprojecting a spherical system (cf. Binney 
& Tremaine 1987). Thus we find a general expression for the phantom spheroid (PS); there is no restriction on the radial 
profile or the axis ratio of the prolate distribution P(|a;|, y/ y 2 + z 2 ), and in fact it is allowed to be oblate. Any PS can also 
lead to a family of PS densities because new PS can be generated by applying a linear operator to an old PS. For example, 
suppose F(x,y,z) has a free parameter j3, then -^F(x,y, z) is also a PS (cf. Appendix B). 

To reformulate the results in the previous section about the "twin bars" in mathematical terms, we define a new set of 
rectangular coordinates (x a , y a , z a ) which relate to the rectangular coordinates (x, y, z) by a rotation around the z-axis (i.e., 
in the xy plane) by an angle a with 

(X a , ya, Za) = (x COS Q + y SU1 OL,—X SU1 Q + y COS Ct, z) . (11) 

The rotation transformation has the property that 

2 2 2 a 2 ■ 2 Oi Ct OL 

x a = x a/2 cos -+y a/2 sm - + 2x a/2 y a /2 sin - cos -, (12) 

2 22 a 22 Q Ct a 

x = x a/2 cos - +y a /2 sin - - 2x a/2 y a /2 sin - cos - . 



In these notations we let P(|x a |, Vr 2 — x' a ) denote any general prolate bar, made by rotating its end-on twin bar P(\x\ , y/r 2 — x 2 ) 
to an angle a from the line of sight in the xy plane. Then we can always design a triaxial object 



T{\x a/2 \, \y a/2 \, \z a/2 \) = P(\x a \, \/r 2 - x 2 ) + P(\x\, yj r 2 - x 2 ) - S{r), (13) 

where S(r) is given by eqs. (^) and ( |Io| ) so that the triaxial object T(\x a / 2 \, \y a / 2 \, \z a / 2 \) has exactly the same surface 
brightness as the prolate bar Pd^cl, yj y a + z a ), i.e., 
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Figure 2. Cross sections of several phantom spheroidal density models (cf. eq. |45|) in terms of x (horizontal axis) vs. \fy' 1 + z 2 (vertical 
axis). The distributions are rotationally symmetric around the Sun-center axis (i.e., the horizontal x-axis). The shaded regions are 
negative density zones of the PS. We show several locations of the observer (solid circle to the right) at Do = 4a, 2a, a, 0.5a (from top to 
bottom) from the object center (solid circle to the lower left corner) and draw a line of sight which intersects the x = plane at a radius 
a. The amount of light adds up to zero along any line from the observer's perspective. 
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/oo p oo 

P(\x a \,y / y 2 +z 2 )dD = / T(\x a/2 \,\y a/2 \,\z a/2 \)dD. (14) 
-oo J — oo 

To verify the mirror symmetries of the distribution T(\x a / 2 \, \y a /2\, \ z a/2\), simply substitute eq. ( |l2| ) to eq. ( |l^ ) and apply the 
transformation x a / 2 — > ~x a / 2 or y a / 2 — > —y a / 2 - In both cases we have \x a \ — > \x\. The mirror symmetry with z a / 2 — plane 
is obvious because T(\x a / 2 \, \y a /2\,\ z a/2\) is an even function of z a / 2 , which comes in only through the r 2 — x 2 a j 2 + y 2 a j 2 + z 2 ^ 
factor. Note that we do not lose generality by choosing the rotation to be around the z-axis since there is no prefered direction 
in the yz plane of the prolate bar P(|a:|, \J y 2 + z 2 ). 

To summarize the main result here, we have shown that any prolate bar has a triaxial counterpart with identical surface 
brightness but half the viewing angle. An immediate question is to what extent we can build triaxial models with the viewing 
angle varying continuously. 



4 A GENERAL SEQUENCE OF TRIAXIAL MODELS WITH IDENTICAL SURFACE BRIGHTNESS 

In the following we will construct a subclass of these phantom spheroids which preserve the mirror symmetry for a continuous 
sequence of fairly general triaxial models. We do not know whether such constructions exist for all triaxial models, but it 
does exist if the triaxial density distribution vo(x,y, z) can be decomposed to a spherical part G(r) and a non-spherical part 
p{r)Q (x,y,z), i.e., 

u (x,y,z) = G(r) +p(r)Q (x,y,z), (15) 
where the only restriction is that Qo(x,y, z) is a quadratic function of rectangular coordinates, i.e., 

Qo = c\\x 2 + c 22 y 2 + C33Z 2 + 2c\ 2 xy + 2c 2 zyz + 2cz\zx. (16) 

There is complete freedom with the functional forms of G(r) and r 2 p(r), which are the the radial profiles of the spherical 
component and the non-spherical component respectively. These models are also triaxial by construction because surfaces of 
constant Qo are ellipsoidal isosurfaces. The symmetry planes of the ellipsoid can be oriented to any direction by changing 
the elements of the symmetric matrix dj where dj = Cj, and the indices i and j are from 1 to 3. The total luminosity of the 
model, Lo, is given by 

/•oo ^ 3 /> oo 

d 3 rvo(x, y, z) = 4m I drr 2 G(r) + — ^^ca I drr A p(r). (17) 
Jo 6 i=1 Jo 

We define 1(1' , b') as the surface brightness integrated along both the (I', b') direction and the opposite (n + I' , —6') direction. 
The central surface brightness, Iq, integrated along the a;-axis (i.e. with x from — oo to oo and y = z = 0) is given by 

/oo 

v (x,Q,Q)dx = / [G(r) + cnr 2 p(r)] dr, r = x, (18) 

i J — oo 

where r and x are interchangeable dummy variables for the integrations. 

For the above triaxial model we can construct a phantom spheroid F(x, y, z) with 

F(x,y, z) = P(\x\, vV + z 2 ) - S(r), P(\x\, y/y» + z») = x 2 p(r), (19) 



where the spherical component S(r) is computed from (cf. eq. [HJ) 

1 f°°dj(w) dw 
n J dw - r 2 



S(r) = -- r dJ ^ *L=, JM = n* 2 P(r)dD, (20) 



so to cancel the P(|a;|, \/ y 2 + z 2 ) component in the projected light distribution. Clearly F(x,y,z) is a spheroidal distribution 
with axial symmetry around the Sun-center axis since it is made by subtracting from a prolate distribution x 2 p(r) a spherical 
distribution S(r) of identical surface brightness. Let Jo be the line of sight integration of either the spherical S(r) or the 
prolate x 2 p(r) distribution along the x-axis, then 



Jo = J(0) = / S(r)dr = / r 2 p{r)dr. (21) 

J — oo J —oo 

The total luminosity of the phantom spheroid, Lps, is given by 

a P OO POO 

Lps = Lp-Ls, Lp = ^-J drr 4 p(r), L s = 4tt J drr 2 S(r). (22) 

Adding a fraction 7 of the above phantom spheroid to our original model vq(x, y, z) we get a new model 

v-y(x,y,z) = is (x,y,z) + F(x,y,z)j, (23) 
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Figure 3. Same as the previous figure but for the new phantom spheroidal densities F\, F 2 and F3, (cf. eq. [t6| ) with the observer at 
D = 4a. 

which should have identical surface brightness map as the old one. Rearranging the terms we find 

v-y(x,V,z) = G 7 (r) +p(r)Qj{x,y,z), (24) 
where 

G 7 = G{r) - jS(r), r = ^x 2 +y 2 + z 2 , (25) 
and 

Q~f{x, y, z) = Qo(x, y, z) + 7s 2 = (7 + cu)x 2 + c 2 2y 2 + c 33 z 2 + 2c 12 xy + 2c 2 zyz + 2c 31 zx. (26) 

So the new model is very similar to the old model. Both are superpositions of a spherical component and a triaxial perturbation; 
the triaxiality is guaranteed by the fact that Q 1 {x,y, z) prescribes ellipsoidal isosurfaces. In general the orientation of the 
object is a function of 7 given by 
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, o 7 + c ll~ c 22 . „ 7 + c ll~ c 33 ,„ C 2 2 - C33 

cot2a xy = , cot2a.« = , cot2ag, z = — , (27) 

2C12 2C13 2C23 

where a IS and a xz are the tilt angles of the object from the the line of sight in the xy plane (i.e., the z = cut) and the xz 
plane (i.e. the y = cut), and a yz is the tilt angle counting from the y-axis in the yz plane (i.e., the x = cut). We make two 
comments here. First the orientation in the yz plane is not affected by the phantom spheroid because the latter is generally 
rotationally symmetric around in the yz plane (cf. eq. |l9| ). Second the directions of the principal axes of the triaxial object 
should generally be eigenvectors of the 3x3 matrix Cij + j8u with i,j = 1, 2, 3. The corresponding eigenvalues are generally 
the three roots of a cubic equation. In practice, it is more convenient to speak of the orientation of the object in terms of the 
angles a xy and a xz , which prescribe the major axes in the xy and xz plane cuts. 

There are, however, limits to 7, the amount of PS that we can add to a model. A physical density model must have a 
positive phase space density and must be stable. This, in principle, set limits on the non-uniqueness if we can build dynamical 
(numerical) models for a sequence of potentials with different 7. In practice stringent limits can already be obtained from the 
minimal requirement that the volume density of a physical model is everywhere positive, i.e., 

u y (x,y,z) = \G(r)+p(r)Q (x,y,z)]+j[x 2 p(r)-S(r)] > 0. (28) 

This generally involves (numerically) searching over a 3-dimensional space (a;, y, z) for the minimum of the density v 7 (x, y, z), 
and finding the range of 7 which brings the minimum above zero; for edge-on models, this can be reduced to a search in 
2-dimensional space (see Appendix A). 

Nevertheless there are many easy-to-use variations of the positivity equation. For example, a positive volume density at 
the object center requires (cf. eq. ^) 

u-,{0, 0, 0) = G 7 (0) = G(0) - 7S(0) > 0, (29) 
while a less interesting condition can come from requiring a positive total luminosity of the object L 7 , given by 

L 7 = Jd 3 rv 1 (x,y,z) = L + -yL PS > 0. (30) 

A stringent set of conditions can be derived by computing the following moments of the density distribution. This is done 
by first multiplying both sides of eq. (^) by a factor r k , where k can be 0, 1, 2 etc.. Then integrate over r along a general 
direction n = (m, 77,2, 773) through the origin, i.e., from r = —00 to r = 00 with (x,y,z) — (nir, 7i2r, 7i3r). We find 

3 poo poo 

v~f(nir, n^r, n 3 r)r k dr = / [G(r) + c nn r 2 p(r)] r k dr + / ^yn,ir 2 p(r) — 7i5(r)l r k dr, (31) 

o J — 00 J — oc 

where 

3 

n\ + nl + n\ = 1, c nn = ajmrij. (32) 

Requiring the k — moment of the density v^{x, y, z) to be positive we have 

3 \ 2 2 



u~f(nir, n 2 r, n 3 r)dr = J + 



Jo > 0, (33) 



where we have substituted in the definitions for Jo and Jo (cf. eqns. ^| and [H]) . Eq. ( |33|) generally yields a necessary upper 
limit on 7, and the limit is most stringent when the line (x, y, z) — (nir, nir, nzr) is the minor axis of the model. For example, 
along the j/-axis with (ni, 712,723) = (0, 1,0) we have 

2ci2 cot 2a xy = 7 + en — C22 < -7^ (34) 

■Jo 

where we express the constraint in terms of the tilt angle et xy (cf. eq. p7[ ) as well as 7. Note we assume Jo > 0, which is 
generally valid. Likewise along the z-axis with (711,712,773) = (0,0, 1) we obtain 

2ci3 cot 2a xz — 7 + en - c 33 < (35) 

Jo 

Likewise a positive k — 2 moment of the density requires 

r 2 ( 3 3 \ 

47r / ^ 7 (7iir, n 2 r,n 3 r)r 2 dr = L + 3 ^ CijUiUj ■- ^ cu \ L P + j(3n(L P - L s ) > 0, (36) 

V i,j=i i=i j 

where we have substituted in the definitions for Lo, Lp and Ls (cf. eq. ^7] and p^) . Eq. (^) generally yields a necessary lower 
limit on 7. For example, along the z-axis with n = (1, 0, 0) we have 
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Figure 4. The run of the phantom densities F (top), F\ (middle) and F2 (bottom) as a function of the radius, where the length is in 
units of the characteristic length a, and the density has units of v c . The run along the z-axis is shown in the left half (the negative half 
of the horizontal axis), and the run along the x-axis is shown in the right half (positive half). Various line styles are for Do = 4a, 2a and 
a with heavier lines for larger Dq. 



2ci2 cot 2a xy + 2ci 3 cot 2a xz = 2j + 2cn - C22 - c 33 > -— ^, (37) 

Lp 

where we assume the prolate component has a positive luminosity, i.e., Lp > 0, which is generally the case. The r.h.s of 
eq. ([37]) can be further simplied in the limit Do — > 00, i.e., the observer is far away from the object. In this case L 1 — > L . 



5 A SEQUENCE OF NEARLY EDGE-ON TRIAXIAL MODELS WITH A GAUSSIAN PROFILE 

The above general results apply to models with any radial profile. To illustrate the model properties effectively, we will use 
the following smooth triaxial models with a Gaussian radial profile with 

vo(x,y,z) = G(r) +p{r)Q (x,y,z), G(r) = v c e 3?*, p( r ) = -|e ~^ , (38) 

where a is the characteristic scale of the model, v c is a characteristic density. These two intrinsic parameters are related to 
the observable 7o (cf. eq. ^) by 
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Io= vo{x, 0, 0)dx = ^ (V2 + 2i) u c a, (39) 

where Io is effectively the integrated light intensity from both the object center (V ,b') = (0,0) and the anti-center (l',b') — 
(7r,0). While we prefer to keep our derivations as general as possible, wherever quantitative numerical calculations are shown 
we set the scaling parameters a and v c to unity, and use the following set of model parameters 
5 7 1 

Cll = C 22 = -g, C33 = -X, Cl2 = -, (40) 

and 

c 3 i = C23 = i, (41) 
and put the object at the distance 

Do = 4a (42) 

unless otherwise mentioned. These parameters describe nearly edge-on models with a characteristic angular size of jy- = 4 ~ 
15°; edge-on models have C31 = C23 = 0. We will vary 7 to show the properties of a sequence of models. The results can be 
generalized qualitatively to all models described by eq. (p4|). 

To find the expression for the phantom spheroid, we substitute eq. (Bq) in eq. (p|). This yields 



./(»-)=/ .!•'■>(,■),//;= ./,„ ^ ( J + ^ - ^ ) , Jo = ^m. m:;:i 



where J(w) is the integrated intensity of the prolate distribution x 2 p(r) along the line of sight with the impact parameter w, 
which becomes Jo when w = 0. Applying eq. ( p"ci| ) we find that S(r) reduces to a simple analytical expression, 

1 f°°dj(w) dw 

nj r dw ^w 2 -r" -■■•■■'\-2 ' u~D~ 2D;: 
and the PS is given by 



S(r) = -zl T , r ^=*Mr)\~ + ^-^), (44) 



F(x, y, z) — x p — S = v c e 



Dl V 2a 2 



(45) 



A few examples of the PS are shown in Fig. ^. They vary as a function of characteristic angular size of the object j^. When 
the observer is well inside the object, 3> 1, the PS density becomes nearly spherical because of the strong dependence on 
r (cf. eq. Wa). The amplitude changes from ~ —77^-^0 to ~ ttti^c when r moves from r — a to r = 1.5a. 



a 



We can also generate a sequence of analytical phantom spheroids by repeatedly applying the derivative operator | gl 
to the original phantom spheroid F(x,y,z) (cf. eq. [45] ), i.e, 



F n ( X ,y, 2) = X * Pn{ r)-S n (r), ff n(r ) s (46) 

for n — 1,2,3, where x 2 p n (r) and S n {r) form a new pair of a prolate bar and a spherical bulge with matching surface 
brightness. These new analytical PS densities are given in the Appendix B. The contours for these new PS densities are shown 
in Fig. |^, and radial profiles in Fig. ^. As Do decreases, the amplitude of the radial oscillation of the PS increases, but at the 
origin we have, 

F(0, 0, 0) = -0.5v c , Fi (0,0,0) = F 2 (0,0,0) = 0, (47) 
independent of Do- 

Adding the phantom spheroidal density F(x,y,z) to our 7 = model (cf. eq. we obtain the general expression for 
the model density 



f~f(x,y,z) = u c e 2^<l + e 



ly _ 7 _ 7r 4 37r 2 
1 2 2 a 2 J> 2 2D 2 

Q 7 = (7 + cu)x 2 + c 2 2y 2 + c 33 z 2 + 2c 12 xy + 2c 2 -j,yz + 2c 31 zx. 



(48) 



The object-observer distance Do can take on any positive value, and the parameters 7 and Cij can describe the most general 
orientation of the three symmetry planes of the model (cf. eq. ^). These models have also the nice property that the PS 
density falls off steeper than the spherical term G(r), meaning that the models are always nearly spherical with a positive, 
Gaussian profile at large radii. 

The free parameter 7 is constrained by the requirement of a positive volume density to the range 

Min(7) < 7 < Max( 7 ), Min( 7 ) ~ -0.5, Max(7) ~ 0.6, (49) 
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Figure 5. Cross sections of a triaxial nearly edge-on model (cf. eq. ^Cj) with the major axis at a xy = 50° from the Sun-center line in 
the xy plane (top) and at a IZ = 23° in the xz plane (bottom). 




Figure 6. Similar to the previous diagram, but with the major axis at a xy = 25° from the Sun-center line in the xy plane (top) and at 
a IZ = 7° in the xz plane (bottom). This model appears identical to the previous model in terms of the surface brightness distribution. 
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1.2 




radius 



Figure 7. The run of the volume density along the x (thickest line), y and z (thinnest line) axes for the a xy = 50° model (upper curves) 
and the a xy = 25° model (lower curves). The curves are blown up at large radius by a factor of 10 and 1000, and they converge to one 
curve because the density at large radius is dominated by the spherical Gaussian component (cf. eq. |48[), 



where the lower and upper limits are found by computing numerically the density on a spatial grid in the (x, y, z) and 
examining the positivity (cf. eq. ^ with the parameters in eq. ffo] ). They translates to 



Min(a xy ) < « xy < Max(a xy ), Min(cv xz ) < « xz < Max(cv xz ), 
where 



Min(cv xy ) 
Max(a xy ) 



1 

- arctan 
2 



2cia 



7T 1 



arctan 



Max(7) + en — c 2 2 
2c 12 



25°, 



C22 — cn — Min(7) 



65°, 



(50) 



(51) 



Min(a xz ) 
Max(a xz ) 



■ arctan 



2 

7T 1 



2c 3 i 



■ arctan 



Max(7) + cn - c 33 
2c 3 i 



7°, 



C33 - cn - Min(7) 



65°. 



Necessary (but not sufficient) limits can also be derived analytically using, for example, conditions eq. (|34[), and (^ 
substitution of eqs. (|43|) and (J39[) for our Gaussian triaxial models we find 



> 



> 



2C12 



1 /2c 12 Jo\ 1 

— arctan I — I = — arctan 

2 \ Io J 2 

1 . /2c 3 iJo\ 1 . 

— arctan = — arctan , _ 

2 V Io J 2 V2V2 + C11 



2^2 + cii 

2C31 



10°, 
3°. 



And likewise eq. (p7|) requires 



7 ( 1 -S) + ( 2 ^ + ^ Cll] ! ' 



(52) 



Upon 



(53) 



(54) 
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For an observer at Do = 4a this reduces to a lower limit on the amount of phantom spheroid with 7 > —1.75, or to an upper 
limit on the model major axis angles with 

a xy = tJ- — - arctan [ — — ) < 80°, (55) 

2 2 \c22-cn-7j 

7T 1 / 2c 3 i , 

a xz — arctan < 85 . 



2 2 \C33 - en - 7 / 

Fig.0 and ^ illustrate two Gaussian triaxial models, where the parameter 7 is fixed by assigning the tilt angle a xy 
(cf. eq. [27p in the xy plane to 25° or 50°. Figs. jlo| ) and ( |l"l| ) show the identical surface brightness distribution of the two 
models, independent of 7. Unlike the cigar-shaped bars in Fig. [j], these two triaxial models are smooth with a positive density 
everywhere and an overall Gaussian radial profile except for a small bump near (cf. Fig. ^). 

The two models differ only in the even part of the volume density as shown by the cross-sections at, say, x = 2a/5 (cf. 
Fig. [s]). Subtracting one model from another we get back the phantom spheroid, which is rotationally symmetric around the 
Sun-object center axis. 

The odd part of the volume density is identical for both models, and is shown in Fig. (|^). Mathematically, the odd part 
of the density is defined by 

v odd, y ^ [ Vl ( x ,y,z) - v-y(x,-y,z)] = 2y (ci2X + C232) p(r), (56) 
and 

v odd,z _ [ v ^{x,y,z) - v^(x,y,-z)] = 2z (ci3X + c 23 y)p(r). (57) 
Clearly both u odd ' v and 

u odd,z are mc jependent of 7, i.e., the amount of PS in the model. 
The odd part of the density shows up as the asymmetry in the surface brightness. The left-to-right asymmetry is a 
line-of-sight integration of v° dd ' v with 

1 f°° 

i[l(l',b')-I(-l',b')] = / v° dd *dD, (58) 
" J — 00 

and the up-to-down asymmetry is an integration of u odd ' z along the line of sight with 
1 f 00 

±[l(l',b')-I(l',-b')]= u odd ' z dD. (59) 

" J —00 

More specific for the Gaussian triaxial models here, the asymmetry along the b 1 — cut is given by 

I{l',0)-I(-l',0) . , . Ci 2 

= I0C12W (w) sin2i , C12 = — 7= , (60) 

2 2V2 + CH 

and the asymmetry along the I' = cut is given by 

7(0, bQ -1(0,-0') . . . , . cia 

^ = 7oCi3W / (w)sm26 , ci 3 = — -= (61) 

2 2v2 + cii 

where C12 and C13 are the combinations of the intrinsic shape parameters of the model, and 

W(w) = e- ! £ (^-l), (62) 



is a function of the impact parameter w. Effectively W(w) is a rescaled asymmetry distribution after removing the periodical 
anti-symmetric term sin 21' or sin 2b' . Note that the asymmetry distribution has a reversal of sign around the impact parameter 

w = a/V2. (63) 

This means we can derive the scale length of the model a directly from the observable asymmetry map. For example, for the 
models in Fig. (|ll]), the object is at a distance Do = 4a, so the reversal of the asymmetry, i.e., W[w) = 0, happens at 

\l'\ =arcsin-^ ~ 10°, (64) 
Do 

in the longitudinal cut, or at \b'\ ~ 10° in the latitudinal cut. While fitting the asymmetry map will tell us unambiguously 
the scale length a and the shape parameter e of the triaxial object, it does not directly constrain the major-axis angle a, nor 
does it constrain the amount of phantom spheroidal density. 
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Figure 8. The cross-sections of the volume density at x = 2a/5 for the a xy = 25° and a xy = 50° models (solid and dashed contours in 
the upper panel). The maxima of the contours are where the the long-axis of the bar intersects with the x = 2a/5 plane. They are offset 
from the y = z = point because the long-axis is pointed away from the line of sight. Subtracting one from another we get the phantom 
spheroidal distribution (lower panel), which has rotational symmetry around the a;-axis. 

6 PHANTOM SPHEROID VS. THE GALACTIC BAR 

This new form of non-uniqueness of nearby objects might apply to the Galactic bar as well. As mentioned in the Introduction, 
the COBE/DIRBE surface brightness distribution is the main source of information about the volume density of the Galactic 
bar. Imaging the effect of adding any amount of the phantom spheroids here ■yF(x,y, z) to the COBE bar. It does not matter 
whether we use the cigar-shaped one in Fig. |l| or the smoother ones shown in Fig. ^, as long as the PS is exactly centered 
at the Galactic center, i.e., we set Dq of these PS densities to the Galactocentric distance (say 8 kpc). The result is a new 
volume density of the Galaxy with a generally twisted shape in the central part, but there should be no trace of alteration in 
the the surface brightness maps (the symmetric maps or the asymmetry maps). This is purely because of the construction of 
the PS (cf. eq.|). So the perspective effect in the COBE/DIRBE maps is bypassed completely. However, adding any phantom 
spheroid here will modify only the even part of the volume density and will have no influence on the anti-symmetric part 



© 1999 RAS, MNRAS 000, 



Deprojection of light distributions of nearby systems 15 




Figure 9. The cross-sections at x = 2a/5 for the odd part of the model density, which is independent of 7. The upper panels shows 
u-y(x, y, 2) — u-y(x, — y, 2), and lower panel shows v~ f (x,y,z) — u-y(x,y, —z). The distributions are anti-symmetric with respect to the 
y = V = plane or the 2 = b' = plane with both positive (thick contours) and negative contours (thin contours). 

since the PS is always symmetric with respect to the I — plane and the b — plane (cf. eq. |^ and Fig. ^ . So as illustrated 
by Fig. (^) and Fig. (p"o|), the prominent perspective effect in the left-to-right asymmetry map [1(1, b) — I(—l, b)[ /2 or the up- 
to-down asymmetry map [1(1, b) — 1(1, —b)[ /2 should tightly constrain the part of the volume density which is anti-symmetric 
with respect to the I = 0° plane or the b — 0° plane (cf. eqs. |1| and p9[ ). 

It is tempting to suggest that phantom spheroids are responsible for the uncertain major axis angle of the Galactic bar 
from the COBE/DIRBE map (Binney & Gerhard 1996, Binney, Gerhard & Spergel 1997). However, this statement depends 
on whether there exists a phantom spheroidal density which preserves the assumed mirror symmetry of the COBE bar exactly. 
Unless the added PS has a density profile matching that of the COBE bar, it will create a m — 2 twist, like a two-armed spiral 
pattern or an S-shaped warp. On the other hand, part of such distortion may be hidden by the spiral arms in the Galactic disk. 
After all the mirror symmetry is merely an assumption, which is motivated by the fact that external face-on barred galaxies 
have a largely bi-symmetric distribution of light. Here we remark that the sequence of Gaussian triaxial models here capture 
the main features of the Galactic bar and the COBE/DIRBE maps: the models shown in Fig. ^ and [] with the scalelength 
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Figure 10. The identical surface brightness map of the two triaxial models from the observer's perspective, where the horizontal axis 
y/D and the vertical axis z/D measure the angular distance of a line of sight from the y = plane and the 2 = plane. The dashed 
contours in the three quadrants are the images of the contours in the first quadrant. The offset with the solid contours shows the 
asymmetry or the perspective effect of a nearby object. Contours are in intervals of 1 magnitude. 



a — Do/i = 2kpc resemble the Gaussian Galactic bar models of Dwek et al. (1995) and Freudenreich (1998) in terms of aspect 
ratios, radial profile and the offset of the Sun from all three symmetry planes of the bar; the asymmetric patterns in the surface 
brightness cuts in Fig. ([n]) also mimics those seen in the COBE/DIRBE maps. Hence it would not be surprising if there exists 
a phantom spheroid which is qualitatively similar the one prescribed by eq. ( |45| ) and preserves the mirror-symmetry of the 
Galactic bar. Nevertheless a quantitative determination of the parameters of the COBE bar is a complex numerical problem 
involving to the least a reliable treatment of dust extinction (Arendt et al. 1994, Spergel 1997). Note that patchy dust in the 
Galaxy makes a phantom spheroid here project to a non-zero irregular surface brightness map. The effect is strongest when 
the dust is mixed with stars near the center, but even a dust screen near the Sun can cast a faint shadow unless the phantom 
spheroid is confined inside the solar circle (e.g. the cigar-shaped PS in Fig. [l]). The fact that dust extinction is a function of 
wavelength could constrain the phantom spheroids. The range of the bar angle and axis ratio are also constrained by other 
observational constraints of the bar such as microlensing optical depth (Zhao & Mao 1996, Bissantz, Englmaier, Binney, & 
Gerhard 1997) and star counts (Stanek et al. 1994). Star count data, for example, can constrain the line-of-sight distribution 
of the bar, hence can in principle distinguish between a prolate bar and a spherical bulge, and break the degeneracy of the 
PS. A detailed numerical treatment of the Galactic bar is beyond the scope of this paper on the existence of non-uniqueness 
in deprojecting a general nearby object. 



7 PHANTOM SPHEROID VS. KNOWN NON-UNIQUENESS IN EXTERNAL GALAXIES 

The phantom spheroidal density for nearby systems is of a different nature from the kind of non- uniqueness associated with a 
simple shear and/or stretch transformation of an ellipsoid. This applies even in the limit that the object is at infinity because 
these transformations normally change the odd part of the density, while the PS density here is strictly an even function (cf. 
eq.|). 



© 1999 RAS, MNRAS 000, §-|l| 



Deprojection of light distributions of nearby systems 17 




-40 -20 20 40 



b' degrees 



Figure 11. the run of the surface brightness along the i'-axis (top) and along the b'-axis (bottom) for our Gaussian triaxial models. 
The solid curves show symmetric part and the dashed curves show the perspective-induced asymmetry (after blown up by a factor of 
four). The I' and b' are the generalized longitude and latitude in units of degrees for the object. 

Our phantom spheroidal density F(x,y,z) is a function of the object distance Do. In general 

F(x,y,z) = F Do ^ 00 (x,y,z)+Da 2 F (r), (65) 

where Fd ^oo is the PS in the limit of Do — ■> oo, and Fa(r) is (up to a constant) the PS in the limit of Do — ■» 0. For a 
hypothetical observer receding from the object, the distance Do increases, the PS density is adjusted slightly to accommodate 
the changing perspective so to preserve this kind of non-uniqueness all the way to extragalactic distance. The common 
ambiguity of an extragalactic axisymmetric bulge with an end-on bar is a very special case of the kind of non-uniqueness for 
nearby objects. 

The divergence of F(x, y, z) at small Do (cf. eq. [35|) has an interesting implication. In the limit that Do — > some 
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regions of the phantom spheroid will have infinitely positive density, and some regions infinitely negative density because 
■yF(x,y, z) — > jDq 2 Fo(r), and Jb(r) changes sign, for example, at r — ay/3/2 (cf. eq. ^B] or Fig. ^J). So to add any finite 
amount of PS would violate the positivity requirement. Hence only the model with 7 = is allowed, and the major axis 
direction is fixed to the direction where the object appears the brightest in the projected map. In other words the kind of 
non- uniqueness discussed here disappears if the observer is very close to the object center. In general the range for 7 becomes 
very narrow if Do < a, and is insensitive to the distance Do to the object if Do > 2a. 

Unlike its external counterpart with Do — > 00, the PS density for a nearby triaxial body is not massless. The total 
luminosity L 7 of our model changes with 7. For the Gaussian model we have (cf. eqs. ^ and p^ ) 

d 3 rv-,(x, y, z) — L + jL P s, (66) 

3 Cn + C22 + C33 
Z 2 H 



3 3 
Lq = n 2 u c a 



2 

Sir 2 u c a 3 \ ( 7a 2 



D 2 

where Lo is the luminosity in the absence of the PS. The PS has a luminosity Lps proportional to Dq~ 2 . 

The normalization u c and the scale a are fixed by the projected light intensity 7o and the asymmetry map W(w) (cf. 
eqs. [ii] and [so) ) , but the total luminosity changes by a fraction of the order of jp. This is because unlike external systems, 
the observed angular distribution of light in a nearby object does not simply sum up to a unique measurement of its total 
luminosity. We may slightly underestimate/overestimate the intrinsic luminosity by about a few percent, which may not be 
easy to detect in reality. Moving away from the object the observer has a full outside view and hence a better determination 
of its luminosity but at the expense of relaxing the constraint on the orientation of the object from the perspective effect and 
positivity. 

The phantom spheroid affects also the dynamics of the model, however, the effect can be fairly difficult to detect. For the 
models in eq. (E3J), the depth of the potential well at the center expressed in terms of the escape velocity V esc is given by 



±V? SC = w c a 2 G{M/L) 



2 / 

- (Cn + C22 + C33 J 



which decreases linearly with increasing fraction (7) of the superimposed PS, where M/L is the mass-to-light ratio. As a 
result, the potential well of the a = 25° model is slightly shallower than that of the more centrally concentrated a — 50° 
model (cf. Fig. |B| and ^|). However, the difference is only 7% in terms of the maximum escape velocity V esc of the models. 
The differences in terms of the mass weighted average velocity dispersion (estimated from the virial theorem) and the circular 
velocity (estimated directly from the potential) are also at only a few percent level, too small to be measured with certainty. 
However it might still be feasible to distinguish the models in terms of the stellar orbits in them. The projected velocity 
distribution should change with the major axis angle. 

Our phantom spheroidal density is also different from konuses since it is not confined to any cone in the Fourier k-space 
even in the limit Do — > 00. Using eqs. (4), (24) and (25) of Palmer (1994) we can compute the Fourier transform of our PS 
(cf. eq. ||). We find 

D(k) = /d 3 rexpHk- r) [p{r)x 2 - S(r)] = Jb -«,(k) + D,7 2 D (k), (68) 



which is a k-space spheroidal distribution around the Sun-center line, where 

3 / 2\ k 2 a 2 k 2 a 2 

Fd ^co (k) = tt5 Vc a A (cos 2 <9 k - - e ~ , (69) 

#k is the angle of the k-vector with the Sun-center line, and 



D f 7 2 D (k)EE ( f-l^ + ^_ )e-^L PS . (70) 



2 „ 2 i. 4 „ 4 \ t 2„2 



7k 2 a 2 k 4 a 



We recover the luminosity of the PS in the limit k — > 0, F(0) = D Do(0) = Lps (cf. eq. B6|). Fig. |I2| shows that even when 



the object is placed at infinity, D(k) is non-zero everywhere except at k= or #k = cos" 1 (cf. eq. |69|). Since any konus-like 
structure (or its triaxial version) must have a certain "cone of ignorance" around a principal axis of the model outside which 
F(k) = (Rybicki 1986, Gerhard & Binney 1996, Kochanek & Rybicki 1996), the kind of non-uniqueness shown in this paper 
has little to do with konuses. The two sequences of models meet only when the object is a face-on or end-on spheroidal system, 
in which case, the PS here is equivalent to a trivial konus with the entire k-space in the "cone of ignorance" . 
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Figure 12. The phantom density in the Fourier k-space, k x vs. k p = \Jk^ + k'j in units of a 1 . The distributions are rotationally 

symmetric around the fc^-axis. The shaded and unshaded zones correspond to negative and positive F(k) respectively. The observer is 
at Do = oo,4a, 1.25a from top to bottom. 
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8 CONCLUSION 

The non- uniqueness in deprojecting the surface brightness maps of a nearby system can originate from the simple fact that a 
spherical bulge can match an end-on prolate bar or a face-on disk in the line-of-sight integrated light distribution. Subtracting 
the spherical bulge from the prolate end-on bar or the face-on disk forms a phantom spheroid (PS), which has both positive 
and negative density regions. The non-uniqueness due to such PS can be characterized by two numbers: the distance (-Do) of 
the object and the amount (7) of the PS density that we can superimpose. By tailoring the PS, it is possible to preserve the 
triaxial reflection symmetries of the model. The orientation of these symmetry planes with respect to our line of sight changes 
as increasing amount of PS is added up to the point when negative density regions appear in the final model. The limits on the 
major-axis angles are given analytically. The phantom spheroidal density here forms a new class of non-uniqueness of entirely 
different nature from known degeneracy in deprojecting extragalactic objects. It does not preserve the total luminosity of the 
system. 

The author thanks Frank van den Bosch and Tim de Zeeuw for many helpful comments on the presentation, and Ortwin 
Gerhard for careful reading of an early draft. 



APPENDIX A: POSITIVITY AS A CONSTRAINT TO THE MAJOR AXIS ANGLE OF EDGE-ON 
MODELS 

Edge-on models are specified by 

C23 = C31 = a yz — a xz — 0, C12 = e, a = a xy , (Al) 

so that the z = plane, which passes through the observer at (x,y,z) = (Do, 0,0), is a symmetry plane of the model. The 
volume density model of the system is best described in a cylindrical coordinate system (R, Z, ip) centered on the object. This 
coordinate system is related to the rectangular coordinate system by 



x — Rcosip, y = Rsinip, z 



(A2) 



where the Z-axis is the symmetry axis of the edge-on system, and if) = is a plane, passing through the Sun-object line and 
the Z-axis. In this coordinate system we can rewrite the ellipsoidal term Q 1 (cf. eq. ^) as 

d2 

Q 7 = WiZ 2 + (c 22 + ecot 2a) R 2 ] + -— - cos(2^ - 2a). (A3) 
L J sm 2a 

Clearly surfaces of constant Q 1 are ellipsoidal surfaces with mirror symmetry with respect to the ip = a plane, the ip = a + ^ 

plane and the Z — plane. 

The triaxial model density (cf. eq. [m|) can then be written in above notations as 

v-j(R, Z, ip) = n (R, Z) + n 2 (R, Z) cos(2V> - 2a), 



where the second term is a triaxial perturbation with the major axis along ip — a and amplitude 

n 2 (R,Z) = . , 
sin 2a 

and the first term, the axisymmetric part, is given by 
no (R, Z) = E(R,Z)+e cot 2a [R 2 p - 25] , 
where 

E(R, Z) = [G+ (cu - c 22 )S] + {c 22 R 2 + c A3 Z 2 )p. 

Imposing positivity requirements to v-/(R, Z, ip) for all values of the azimuthal angle tp, we have 
n (R,Z) > n 2 (R,Z), 



i.e., 



eR 2 p 



E(R, Z)+e cot 2a \R 2 p - 2S] > 

L J sin 2a 

Using the following relations for sinusoidal functions 



l + r 



cot 2a - 



t = tana, 



sin 2a 2t ' 2t 
the inequality reduces to an upper and a lower bound for the angle a with, 

Max [t_] < tana < Min [t+] , 



(A4) 

(A5) 
(A6) 
(A7) 
(A8) 

(A9) 

(A10) 
(All) 
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where t- < t < t+ is the range bounded by the effectively quadratic inequality for t 
[eR 2 p -eS]t+ eS^ 1 < E(R, Z) 



(A12) 



at a given position (R, Z) on the meridional plane, and the overlapped interval of these ranges is used in eq. (All). For 
example, for models with parameters in eq. ( [Io| ) and Do = 4a, we find 77° > a > 15° to ensure a positive density at 
(R,Z) = (a, a), but a narrower range 63.5° > a > 22.5° is necessary to ensure positivity everywhere in the (R,Z) plane. 

Interestingly we recover the results of eq. (^ij) if we integrate both sides of eq. (A12) along the y-axis. And likewise if we 
multiply both sides of eq. (A12) by a factor x 2 and then integrate along the x-axis, we recover the results of eq. (pS7|). 



APPENDIX B: ANALYTICAL EXPRESSIONS FOR A SERIES OF PHANTOM SPHEROIDS 

New phantom spheroids can be generated from the original PS F(x,y, z) = x 2 p(r) - S(r) (cf. eq. ||) by repeatedly applying 
the derivative operator 



2 S In 



Fi(x,y,z) = dF ^: V '* ] = x 2 Pl (r) - S 1 (r), 



where 
Pi(r) = 



2<9(lna) 



dp{r) , , ( r 2 . 
— = -p(r) 1 J , 



2 2 

r a 



2a 4 



and 

^ / x dFi(x,y,z) 2 / \ , x 
F ^(x,y,z) = _ . ' ' = x 2 p 2 (r) - S 2 (r), 



2<9(lna) 



where 
P2(r) = 

S 2 (r) = 



dpi(r) 
2d (In a) 

<9Si(r) 
29(lna) 



p(r) 
- —a 2 p(r) 



3r 



1 r l 

2 + „4 



2a 2 2a 4 



+ 



D 2 



5r 4 9r 6 
"2^ + Ifi 



(Bl) 

(B2) 
(B3) 

(B4) 

(B5) 
(B6) 
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